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Numerical Determination of Short Period Trojan Orbits 
in the Restricted Th^ee Body Problem 

Edson F. Goodrich 
ABSTRACT 

In the plane restricted three body problem, two classes of periodic 
orbits exist around the equilateral libration points. Rabe developed 
methods for determining members of the class of long period orbits 
with Jupiter and the Sim as principal masses. This report demon- 
strates the use of these methods to examine the class of short period 
orbits. Power series expansions permit the solution of the relevant 
equations of motion on an IBM 7094 computer; closed form recurrence 
formulae allow the computation of the coefficients of the power series. 

Knowledge of the value for the Jacobi constant gives approximate initial 

♦ 

conditions for short period orbits near the libration points. An itera- 
tion scheme improves the initial conditions to periodic conditions. For 
the starting positions used, two types of short period orbits are com- 
puted with an accuracy of twelve significant figures. The Jacobi con- 
stant decreases as the orbit size increases and deviates farther from 
the libration point. The periods of all orbits differ by less than one 
percent from the period of Jupiter’s motion around the Sim. Elliptical 
orbits around the Sun closely represent the short period orbits, and 
the eccentricities of these elliptical orbits approach unity as the devia- 
tions from the libration point increase. A linear stability study 
• > 

indicates instability for all orbits determined. 
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i. INTRODUCTION 


The restricted problem of three bodies involves the study of the motion of 
.ui infinitesimal particle moving under the influence of two finite point masses 
and is one of the classical problems in celestial mechanics. Studies show the 
ex'sumcc of periodic orbits of the immitesimal mass about the two equilateral 
Lagrangian libration points. These periodic orbits fall into two general classi- 
fications consisting of short period and long period orbits. The short period 
orbits have periods on the order of the period of the two principal masses about 
their common center of mass. The periods of long period orbits depend on the 
mass ratio of the finite masses. With the Sun and Jupiter as predominant masses, 
these orbital periods are on the order of twelve times the period of Jupiter’s 
motion around the Sun. 

Most numerical investigations give results of a few significant figures. 

Rabe (1961) showed to a high degree of accuracy that orbits exist for the class 
of long period orbits. To a limited degree of accuracy, Willard (1913) examined 
the short period orbits. This present work is an effort to characterize, to a 
high degree of accuracy, the family of short period orbits. 

Steffensen (1956) developed a numerical integration scheme for solving the 
equations of motion in winch he introduced two auxiliary variables and expanded 
the coordinates in powers of the time. He derived recurrence formulae in closed 
form for the coefficients of the power series. 

For the case considered, the Sun and Jupiter represent the predominant 
masses and a Trojan asteroid represents the infinitesimal mass. By assumption, 
the Sun and Jupiter move about their common center of mass in concentric circles 
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undisturbed by the Trojan's presence and the Trojan remains in the plane of the 
orbits of the Sun and Jupiter. 

The x-y reference system used has its origin at Jupiter and rotates with 
uniform angular velocity with the Sun fixed in the system on the positive x-axis. 
The unit of distance is the Jupiter-Sun distance, the unit of mass is the mass of 
the Sun, and the unit of time is such that the gravitational constant is unity. If 
M, N, and P represent the mass, mean motion, and period of Jupiter respec- 
tively, then N 2 = 1 + M and P * 2 tt/N. For this study, M = .00095478610 from 
Hill's value for the mass ratio of Jupiter and the Sun. 

Modification of one of the auxiliary variables of Steffensen's method by the 
factor M, as suggested by Babe (1961), speeds convergence of the power series. 
Steffensen's method is ideally suited for use with a high speed computer. Ap- 
proximate starting conditions for a periodic orbit of short period used with an 
iterative procedure yield periodic orbits of high precision. This report considers 
orbits around libration point L 5 . Symmetry considerations yield results for orbits 
around the other equilateral libration point L 4 . 

2. DETERMINATION OF APPROXIMATE INITIAL CONDITIONS 

Finding periodic orbits by an iteration procedure requires starting values 
x 0 , y 0 , * 0 , y 0 which result in a return of the particle to the near vicinity of the 
starting position with velocity components of the same order of magnitude and 
algebraic sign as the initial values. To find the initial conditions for the long 
period orbits, Rabe assumed that, at a point A, the first and second derivatives 
of the velocity with respect to time were zero. This starting point on the line 
joining the Sun and L s is such that L s lies between the Sun and the point A as 
shown in Fig. 1. This assumption enables the derivation of a quadratic 
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Figure 1— The rotating coordinate system 
with a typical periodic orbit. 


expression for the initial velocity components in terms of position. The smaller 
of the two roots to this equation proves satisfactory for the approximate long 
period orbit starting conditions. The larger root to the quadratic should yield a 
good approximation for starting the short period orbits, but trial results proved 
negative. The large first and second derivatives for those orbits invalidate the 
quadratic expression. 

For long period orbits, the velocity at point A is nearly perpendicular to the 
line SA connecting the Sun and the point A. This should hold for short period 
orbits which deviate only slightly from the libration point. The final results 
verify this. If P represents the ratio of velocity components at the starting 
point such that p~ y 0 /x 0 at t = 0, then p is approximately . 577 for small orbits* 
Furthermore, for short period orbits, Charlier (1905), shows that the Jacobi 
integral yields C < 3 (1 + M) where C is the Jacobi constant. This inequality 
relates the initial position to the initial velocity, v - (x 0 2 + y 0 2 )*. These two 
conditions permit the determination of good approximate initial velocity com- 
ponents for starting positions along the line SA and near to the libration point. 
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3. ITERATION TO A PERIODIC ORBIT 


An IBM 7094 digital computer performs the computations involved in double 
precision. The approximate starting position components x 0 , y 0 and velocity 

0 

components x„ , y 0 are read into the computer along with the time interval and 
tolerance desired for the numerical integration. The step-by-step integration 
terminates when the trajectory approaches the starting point A and crosses the 
line SA. The point of crossing is found to a high degree of accuracy, the final 
time interval adjusts such that, within the tolerance specified, the final point of 
the orbit coincides with a point on SA. The orbit does not close sincetheinitial 
conditions approximated periodic conditions. We must .now improve the orbit. 

Holding the starting position fixed, we must alter the initial velocity such 
that the orbit is periodic. If the subscript f denotes. the final values of the com- 
ponents, then y f = y(x f ) since the end position is forced to lie on SA. From the 
Jacobi integral, C = C(x f ,y f ,x f ,y f ), and for points along SA greater than one 
unit from the Sun, C is monotone increasing. Therefore, if x f and y f are equal 
tox 0 andy 0 respectively, x f = x ? andy f -y 0 since C remains constant. A 
differential correction procedure proves satisfactory for improving the initial 
conditions for orbits with small values of A., where X. is the distance from A to 
the libratien point. This method fails for values of A. greater than approximately 
.50. 

A brute force method used in this difficult region actually works for any 
value for X.. With this method, the velocity components for a given \ axe specified 
in terms erf p and v. A small arbitrary change mads in p while holding v fixed 
permits the computation of the change in the quantity Ax = x f - x 0 with respect 
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to (). Linear interpolation gives a new value for p which forces Ax toward zero. 

Ax decreases by several orders of magnitude rather than become zero since the 
system is nonlinear. Repetition of the procedure reduces Ax to the required , 
tolerance. Now introduce a small arbitrary change in y. This causes A* to 
become large again and p must again be adjusted until Ax meets the tolerance. 
These two sets of starting conditions with Ax arbitrarily smti.ll permit the com- 
putation of a new value for v which forces Ay = y f - y 0 toward, zero. Again, the 
reduction is by several orders of magnitude. Iteration of this process' forces 
Ax and Ay to meet the tolerance. This procedure worksJ>egk if the arbitrary 
changes in-p and v. cause Ax and Ay to change sign and decrease in absolute value* 

Computing the Jacobi constant and noting any resulting change provides a 
continuous check cm the computational accuracy , Variations in C indicate corres- 
ponding errors in the accuracy of the results. The converse is not true,, however 
since C can re main "constant though errors exist in the results. Solving the 
equations of motion using a high order Runge-' .utta integration scheme provides 
a final check on the results. Examination of the tinal periodic orbits in this 
manner shows all orbits correct to at least twelve significant figures. * 

A. FINDINGS ON THE SHORT PERIOD ORBITS - 

Periodic orbits were found for values of k beginning with K = . 02 and in- 
creasing in equal intervals of . 02 up to .50. When plotted in the rotating 
frame of reference, the first orbits determined appear to be elliptical in shape 

u 1 

" as shown in Fig. 2. The shape changes as A. increases beyond a value of about 
. 24 where the side nearest the Sun flattens out and eventually bends inward 
toward L s and away from the Sun.' As A. increases further, the orbits dip closer 
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Figure 2— Short period orbits for 
X= .12, .24, .36, .48. 


and closer to the x-axis in the vicinity of the libration point L 3 as illustrated in 
Fig. 2. Finally, as seen in Fig. 3, the orbits ci^ss the x-axis. 



Figure 3— Shc. t period orbit for X = .50 
and the branch point *bit. 


An attempt to iterate initial conditions for X = . 52 to a periodic orbit failed. 
Further examination revealed that two different periodic orbits exist for X ~ . 514, 
no periodic orbit of this type exists for X = . 515 and indeed two orbits exist for 
value of X below some branch point value where the orbits are identically the same 
orbit. This branch point orbit, shown in Fig. occurs at X = .514 325 370. 
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As decreases L*orr. the branch point, the second type of orbit dips further 
below the x-axis and passes closer and closer to the Sun. Passing near the Sun 
results in large coefficients of the power series expansions; and for passage too 
close, the coefficients exceed the computer's capacity. A reduction of the time 
interval decreases the number of terms required in the expansion; however, the 
powers of the time interval become too small for the machine. As a result, we 
show the second type of orbit for values of A. oily down to X = . 36. Figure 4 
represents several, of these orbits. 



Call the smaller of the two orbits for a given X Type I, and call the larger 
Type n. Figure 5 identifies the two orbits for X = .40. All Type II orbits cross 
the x-axis. Type I orbits for X ** . 48 and less do not cross the x-axis, but for 
A ~ .50 and greater they do cross. For X = .496 690 858, the orbit just touches 
the x-axis in the neighborhood of x = 1.85. 
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5. STARTING CONDITIONS AND ORBIT PERIODS 

Table I and Table III give initial position and velocity components for Type I 
and Type n orbits respectively. The x 0 and y 0 coordinates depend linearly on 
a.. Both c omp onents of velocity increase for Type I orbits as X increases to the 
branch point. For Type H orbits, as X decreases from the branch point, x 0 
decreases but y increases. The branch point values for quantities in these 
tables are Xo = .242 837 315, y 0 * 1.311 444 240, Xq = . 749 352 314, 
y 0 = .829 243 523. 

Tables II and IV list the values fear the period T and the Jacobi constant C 
as well as P and v for Type I and Type H orbits respectively. Branch point 
values are p * 1.106 613 681, v •« 1. 117 664 400, T « 6.284 760 760, 

C = 2.367 857 918. Figure 6 shows a plot of p vs X for short period orbits. 
p increases nearly linearly with X for Type I orbits until X = .40. Figure 7 is a 
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Figure 6 — Variation of with \ for 
periodic orbits. 



similar plot of v vs X. Again the behavior is nearly linear until X * .40. 
v roaches a maximum (at least a local maximum) of just over 3 . 20 in the 
neighborhood of X = .44 for the Type n orbits. 

The period and Jacobi constant decrease with increasing X for Type I orbits. 
They continue to decrease as k decreases from the point for Type n orbits. 
Notice that in the vicinity of X - .42 there exists an orbit for which the period 
equals the period of Jupiter, P = 6. 280 187 905. 
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In a nonrotaling frame of reference these short period Trojans move In 
perturbed elliptical two-body orbits with the Sun as the principal body and with 
Jupiter supplying the disturbing force in a periodic manner. An ellipse will, 
therefore, approximate the Trojan motion. Using the initial position and velocity 
in the fixed reference system, we compute the semi-major axis, a, and the ec- 
centricity, e, as the approximating elliptical elements for each orbit. Table V 
and Table VI list a and e for Type I mid Type H orbits respectively. For the 
branch point orbit, a = 1.001 076 088 and e - .730 135 639. . 

6. COMPARISON WITH LONG PERIOD ORBITS 

Notice the slight variations of the periods over the entire range of orbits 
computed. The period deviates less than one-third of one percent from the 
period of Jupiter. Contrast this with the long period orbits which range from 
about six times the orbital period of Jupiter to an orbit of infinite period. 

For Type I short period orbits, the value of the Jacobi constant approaches 
the quantity 3(1+M) from below as X approaches zero. For the long period orbits, 
C tends toward the same value from above as X goes to zero. A plot of the long 
period orbit with the short period orbit for X = .02 in Fig. 8 shows the relative 
sizes of the two orbits. 


Figure 8— Short and long period orbits 
for x = .02. 


0 0.1 0.2 0.3 0.4 0.S 0.0 0.7 0.3 0.0 t.O 
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As * increases, the long period orbits take on a horseshoe shape and are 
<r*ite different from the short period orbits which only gradually bend around the 
>‘»i: . it appears that, as 4 decreases the Type n orbits will take on the horseshoe 
shape, though they will be much thicker than the long period orbits and not sym- 
metrical with respect to the x-axis. 

i or the short period Di bits, the initial velocities exceed those for the corres- 
ponding long period orbits. This accounts for the shorter periods and lower values 
of the Jacobi constant for the short period orbits. 

7. LINEAR STABILITY 

To examine the stability of the short period orbits, we use Hill’s first order 
equation 


- »• 

as treated by Message (1959) and Ra.be (1961>. In Hits equation, rj represents 

the transversal displacement of a Trojan deviating from a periodic reference 

* 

orbit in some disturbed trajectory. If m - then u = m defines the new 
independent variable. Knowing a set of special values for points along the orbit, 
we expand the function 0(u) in a Fourier series of the form 


0(u) 



1 
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Assume the solution of Hill's equation has the form 


T? 


L 


■n k e i(k+c)u 


Substitution of the expressions for 6(u) and r, into Hill's equation produces a set 
of identities from which successive approximations for c can be made. 

From the assumed form of the solution, notice that if c is real, stability 
exists since v must be periodic and the transversal displacement is bounded. 
However, for c complex, the orbit is unstable since 17 then increases without 
bound either when time increases or decreases to infinity. 

The orbits appear unstable since three approximations to c all prove com- 
plex. In all cases, the second and third estimations for c agree to two significant 
figures indicating probably convergence of the sequence. This indicated insta- 
bility for the first order study, especially for the oxbits of small values of a, 
leads to contradictory conclusions. For infinitesimal orbits, the actual motion 
of a non-periodic Trojan will never deviate substantially from the equilateral 
Iteration point. Therefore, at least for small orbits, the motion should be 
stable. It appears that for the short period orbits the oscillation of r/ does not 
remain infinitestimally small, and large values of 17 destroy the first order 
accuracy. A higher order study should clarify this stability problem. 

8 . CONCLUSION 

Several members of this family of short period Trojan orbits closely re- 
semble those found by Willard except for larger values of K; however, he himself 
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questions the results for his large orbits. In this present work, sixteen signifi- 
cant digits were carried throughout the computations and a minimum of twelve 
figures of accuracy were obtained. The Jacobi constant, computed at each step 
of the integration agrees to at least eleven significant places. 

As a. decreases from the branch point for Type II orbits, the general trend 
of the family is well indicated. The orbits pass closer and closer to the Sun 
while bending further around it. SHU, it is difficult to visualize a natural ending 
to the family. Collision with the Sun is the most probable. This question is left 
for further study. 
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Table I 


Initial Conditions for Type I Short Period Orbits 










Table n 


Values of , v , T, C for Type I Orbits 


.* 

/ 

V 

T 

C 

.02 

.590 844 686 

.039 738 614 

6.300 603 307 

3.002 469 808 

.04 

.604 656 383 

.079 113 360 

6.300 560 685 

3.001 283 555 

.06 

.618 387 438 

.118 160 867 

6.300 489 448 

2.999 297 572 

.06 

.632 054 922 

.156 917 640 

6.300 389 288 

2.996 498 052 

.10 

.645 677 281 

.195 420 360 

6.300 259 756 

2.992 865 000 

.12 

.659 274 632. 

.233 706 228 

6.300 100 256 

2.988 371 681 

.14 

.672 869 115 

.271 813 357 

6.299 910 032 

2.982 983 933 

.16 

.686 485 339 

.309 781 235 

6.299 688 156 

2.976 659 259 

.18 

.700 150 929 

.347 651 280 

6.299 433 511 

2.969 345 680 

.20 

.713 897 237 

.385 467 535 

6.299 144 770 

2.960 980 227 

.22 

.727 760 256 

.423 277 530 

6.298 820 369 

2.951 486 985 

.24 

.741 781 825 

.461 133 395 

6.298 458 466 

2.940 774 502 

.26 

.756 Oil 259 

.499 093 326 

6.298 056 899 

2.928 732 318 

.28 

.770 507 569 

.537 223 555 

6.297 613 110 

2.915 226 242 

.30 

.785 342 605 

.575 601 077 

6.297 124 059 

2.900 091 769 

.32 

.800 605 580 

.614 317 529 

6.296 586 078 

2.883 124 701 

.34 

.816 409 835 

.653 484 894 

6.295 994 677 

2.864 067 354 

.36 

.832 903 333 

.693 244 236 

6.295 344 227 

2.842 587 507 

.38 

.850 285 716 

.733 779 716 

6.294 627 459* 

2.818 244 797 

.40 

.868 837 720 

.775 342 438 

6.293 834 611 

2.790 433 880 

.42 

.888 975 849 

.818 294 210 

6.292 951 849 

2.758 280 861 

.44 

.911 364 883 

.863 196 419 

6.291 958 085 

2.720 434 612 

.46 

.937 185 623 

.911 018 660 

6.290 817 485 

2.674 581 000 

.48 

.968 937 726 ! 

.963 755 318 

6.289 457 405 

2.616 020 220 

.50 

1.014 247 299 



1.027 301 748 

6.287 665 797 

2.531 098 829 
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Table m 


Initial Conditions for Type II Short Period Orbits 


X 

r~ ■ 1 

x o 

y 0 


x o 



yo 


.51 

.245 

000 

000 

1.307 

698 

360 

.750 

773 

138 

.875 

606 

527 

.50 

.250 

000 

000 

1.299 

038 

106 

.744 

193 

944 

.909 

782 

831 

.49 

.255 

000 

000 

1.290 

377 

852 

.735 

677 

484 

.930 

854 

592 

.48 

.260 

000 

000 

1.281 

717 

598 

.726 

306 

297 

.946 

762 

135 

.47 

.265 

000 

000 

1.273 

057 

344 

.718 

385 

031 

.959 

686 

855 

.46 

.370 

000 

000 

1.264 

397 

090 

.706 

053 

775 

.970 

609 

262 

.45 

.275 

000 

000 

1.255 

736 

835 

.695 

391 

653 

.980 

070 

591 

.44 

.280 

000 

000 

1.247 

076 

581 

.684 

449 

000 

.988 

407 

217 

.43 

.285 

000 

000 

1.238 

416 

327 

.673 

260 

456 

.995 

845 

229 

.42 

.290 

000 

000 

1.229 

756 

073 

.661 

851 

234 

1.002 

545 

323 

.41 

.295 

000 

000 

1.221 

095 

819 

.650 

240 

472 

1.008 

626 

625 

.40 

.300 

000 

000 

1.212 

435 

565 

.638 

443 

171 

1.014 

180 

390 

.39 

.305 

000 

000 

1.203 

775 

311 

.626 

471 

402 

1.019 

278 

400 

.38 

.310 

000 

000 

1.195 

115 

057 

.614 

335 

072 

1.023 

978 

354 

.37 

.315 

000 

000 

1.186 

454 

803 

.602 

042 

457 

1.028 

327 

485 

.36 

.320 

000 

000 

1.177 

794 

549 

.589 

600 

559 

1.032 

365 

065 
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Table IV 

Values of p , v , T, C for Type n Orbits 


k 

p 

V 

T 

C 

.51 

1.166 

273 

116 

1.153 

406 

735 

6.283 

419 

569 

2.277 

381 

564 

.50 

1.222 

507 

706 

1.175 

384 

799 

6.282 

467 

090 

2.204 

918 

285 

.49 

1.265 

302 

544 

1.186 

470 

240 

6.281 

900 

834 

2.157 

773 

902 

.48 

1.303 

530 

121 

1.193 

264 

169 

6.281 

486 

197 

2.120 

965 

157 

.47 

1.339 

624 

384 

1.197 

583 

555 

6.281 

158 

537 

2.090 

320 

805 

.46 

1.374 

695 

946 

1.200 

247 

589 

6.280 

888 

762 

2.063 

941 

724 

.45 

1.409 

379 

285 

1.201 

710 

412 

6.280 

660 

813 

2.040 

765 

689 

.44 

1.444 

091 

842 

1.202 

255 

905 

6.280 

464 

722 

2.020 

123 

410 

.43 

1.479 

138 

155 

1.202 

076 

271 

6.280 

293 

806 

2.001 

558 

590 

.42 

1.514 

759 

316 

1.201 

309 

361 

6.280 

143 

329 

1.984 

742 

094 

.41 

1.551 

159 

407 

1.200 

058 

473 

6.280 

009 

797 

1.969 

426 

214 

.40 

1,558 

521 

009 

1.198 

403 

750 

6.279 

890 

547 

1.955 

418 

229 

.39 

„ 1.627 

015 

052 

1.196 

400 

158 

6.279 

783 

497 

1.942 

564 

151 

.38 

1.666 

807 

577 

1.194 

126 

982 

6.279 

686 

988 

1.930 

738 

219 

.37 

1.708 

064 

726 

1.191 

600 

830 

6.279 

599 

668 

1.919 

835 

830 

.36 

1.750 

956 

726 

1.188 

867 

717 

6.279 

520 

427 

1.909 

768 

630 
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Table V 


/pproximate Elliptical Elements for Type I Orbits 

a | e 


.02 

1.001 

088 

026 i 

.018 

895 

816 

• 04 JW 

1.001 

215 

557 i 

.038 

773 

078 

.06 ° 

1.001 

338 

042 

.058 

705 

611 

.08 

1.001 

455 

229 ; 

.078 

723 

303 

.10 

1.001 

566 

852 i 

.098 

857 

501 

.12 

1.001 

672 

630 ! 

.119 

141 

244 

.14 

1.001 

772 

255 i 

.139 

609 

582 

.16 

1.001 

865 

389 

.160 

299 

685 

.18 

1.001 

951 

656 i 

.181 

252 

776 

.20 S 

1.002 

030 

633 i 

.202 

511 

962 

.22 1 

1.002 

101 

842 i 

.224 

125 

930 

.24 

1.002 

164 

741 ! 

.24 o 

148 

665 

.26 ! 

1.002 

218 

706 1 

.268 

641 

255 

.28 

1.002 

263 

019 

.291 

673 

953 

.30 i 

1.002 

296 

844 ! 

.315 

329 

050 

.32 

1.002 

319 

195 

.339 

704 

979 

.34 

1.002 

328 

893 i 

.364 

922 

354 

.36 i 

1.002 

324 

502 

.391 

133 

226 

.38 

; 1.002 

304 

221 

.418 

535 

918 

.40 ! 

1.002 

265 

708 

.447 

400 

242 

.42 

1.002 

205 

770 

.478 

113 

729 

.44 

1.002 

119 

722 

.511 

275 

355 

.46 

1.001 

999 

921 

.547 

915 

260 

.48 

1.001 

831 

487 

.590 

143 

082 

.50 

1.001 

572 

690 

.644 

166 

981 
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Table VI 


Approximate Elliptical Elements for Type II Orbits 


a 


.51 

1.000 817 

.50 

1.000 623 

.49 

1.000 504 

.43 

1.000 415 

.47 

1.000 343 

.46 

1.000 283 

.45 

1.000 233 

.44 

1.000 189 

.43 

1.000 150 

.42 

1.000 116 

.41 

1.000 086 

.40 

1.000 059 

.39 

1.000 036 

.38 

1.000 014 

.37 

.999 995 

.36 

.999 978 


e 


788 

.769 986 921 

785 

.798 637 191 

385 

.815 877 208 

094 

.828 623 875 

471 

.838 781 100 

851 

.847 205 780 

066 

.854 371 430 

133 

.860 572 473 

706 

j .866 OOo 795 

821 

.870 814 967 

763 

i .875 101 074 

980 

i .878 944 708 

038 

.882 408 323 

588 

i .885 541 970 

343 

1 .888 386 466 

066 

I .890 975 601 
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Figure 


Caption 


1 The rotating coordinate system with a typical periodic orbit. 

2 Short period orbits for • ~ .12, .21. .30, .-is. 

3 Short period orbit for \ = .50 and the branch point orbit. 

-1 Tyne II short period orbits for - - . 4s, .42, .30. 

5 Type I awl Type 11 short period orbits for • = . io. 

6 Variation of with \ for periodic orbits. 

7 Variation of : with \ for periodic orbits. 

8 Short and long period orbits for \ = . 02. 
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